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abstract 

Using large N, 4d anomaly induced one-loop effective action for conformally invariant 
matter (typical GUT multiplet) we study the possibility to induce the primordial spheri- 
cally symmetric wormholes at the early Universe. The corresponding effective equations 
are obtained in two different coordinate frames. The numerical investigation of these 
equations is done for matter content corresponding to A/" = 4 SU{N) super Yang-Mills 
theory. For some choice of initial conditions, the induced wormhole solution with increas- 
ing throat radius and increasing red-shift function is found. 
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The wormholes are puzzling topological objects (kind of handles of topological origin) 
which attract much attention in General Relativity for years. It is not strange as they 
may be considered as bridges joining two different Universes or two separate regions of 
the same Universe (for an introduction, see refs.[|I|). There are many speculations related 
with hypotetical effects which may be expected near wormholes. 

Moreover, it is known that wormholes usually cannot occur as classical solutions of 
gravity (with matter) due to violations of energy conditions Nevertheless, one can 
expect that primordial wormholes may present at the very early Universe where quantum 
effects play an essential role (or yet in fundamental M-theory, see,for example, 0). 

Indeed as it was shown in ref. (see also |Q) using quantum stress tensor for confor- 
mal scalar on spherically symmetric space and in ref . using one- loop effective action in 
large N and s-wave approximation for minimal scalar there may exist semiclassical quan- 
tum solution corresponding to a Lorentzian wormhole connecting two asymptotically flat 
regions of the Universe. The corresponding spherically symmetric wormhole solution has 
been found numerically in both approximations Q as well as analytically [Q. That 
shows the principal possibility of inducing primordial wormholes at the early Universe (in 
its quantum regime). However above discussion has been limited strictly to scalar matter. 
And what happens if other types of matter are included? Hence, the very natural question 
is: Can primordial wormholes be induced from GUTs at the early Universe? 

In the present Letter we try to answer this question for spherically symmetric worm- 
holes. We use 4d conformal anomaly induced effective action in one-loop and large N 
approximation so that one can neglect quantum gravitational contributions. The effec- 
tive equations for an arbitrary massless GUT containing conformal scalars, spinors and 
vectors are explicitly obtained for two forms of spherically symmetric background. Nu- 
merical solution of these equations is presented for A/" = 4 super Yang-Mills theory. It 
is observed that it depends on the choice of the initial conditions. The initial conditions 
admitting inducing of primordial wormholes are also presented. 

We first derive the effective action for conformally invariant matter (for a general 
review of effective action in curved space, see0). Let us start from Einstein gravity with 
Nq conformal scalars Xi? vectors and N1/2 Dirac spinors ipi 

1 ^ ^ I 1 -^^0 

+ - -A ^ + ^ ^^^^^ • 

" ^ ^J=l k=l J 

The above matter content is typical for asymptotically free GUT at high energies |^ 
as interaction terms and masses at strong curvature are negligible due to asymptotic 
freedom. For finite GUTs or asymptotically non -free GUTs one should not include 
masses if consider only conformally invariant theories and interaction plays no role as 
anyway we consider purely gravitational background. 

The convenient choice for the spherically symmetric space-time is the following: 

ds^ = f{<P)\r\<P)9f.udx^'dx' + rldQ\ . (2) 
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where /i, = 0, 1, g^i, and /(0) depend only from and Tq is non-essential constant. 

Let us start the calculation of matter effective action on the background (0) closely 
following to ref. 0(see also 0). In the calculation of effective action, we present effective 
action as : T = Ti^d + T[l,g^^}] where Ti^d = ^[f,gl^J] - ^[^,9^^)] is conformal anomaly 
induced action which is quite well-known , gj^J is metric (0) without multiplier in front 
of it, i.e., gj^J corresponds to 



g^ydx^dx" + rldVt 

The conformal anomaly for above matter is well-known 
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where b 



(jVo+6Ari/2+12Afi) 



F + -nRj+b'RGB- 

(Wo+llA^i/2+62Afi) ,„ 



b"nR 



(3) 



(4) 



but in principle, b" may be 



120{47rp ' 360(47r)2 

changed by the finite renormalization of local counterterm in gravitational effective action, 
F is the square of Weyl tensor, Rgb is Gauss-Bonnet invariant. 

Conformal anomaly induced effective action Find may be written as follows M: 



W 



b I d'^x^/^Fa + b' / d^x^/^[ 

cr + (Rgb 



a 



-^Rn + '^{V^R)V, 



2 



1 



12 



-77;{b" + 7;{b + b': 



d x\ 



-9 



{R - GDa - 6(V(t)(V(t)}^ - i?^ 



(5) 



where a = |ln/(0), and cr-independent terms are dropped. All 4-dimensional quantities 
(curvatures, covariant derivatives) in Eq.(^ should be calculated on the metric (|^). The 
calculation of r[l, ^fj^^^^] is done in ref. in leading order as follows: 



nh9^'^ 



d^Xy 



bF + b'RcB + 



(6) 



where is mass- dimensional constant parameter, all the quantities are calculated on the 
background (^. The condition of application of above expansion is |-R| < R^ (curvature 
is nearly constant), in this case we may limit by only few first terms. 

Let us solve the equations of motion obtained from the above effective Lagrangians 
S* -|- r. In the following, we use g^y and cr as a set of independent variables and we write 
g^iu as g^y if there is no confusion. 

Lmd {W in Eq.(|^)) is rewritten after the reduction to 2 dimensions as follows: 



J- ind 

An 



^-^Id^x 



+b'rlj d^x^^a (^2n2 + 4i?(2)M-v^V, - ^{R^^^ + Rn)^ 



2 

3 
1 

12 



cr 



+ 3 + 



d Xy 



-9 



X 



{r^'^ + R, 



QUa - 6V^crV„or 



R^'^ + R^Y] 



(7) 
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Here Rq 



^ is scalar curvature of S"^ with the unit radius. 



'^0 



Let us derive the effective equations of motion. In the following, we work in the 



conformal gauge : g±^ 



ip2p 



g±± = after considering the variation of the effective 



action F + with respect to g^i, and a. Note that the tensor g^^^, under consideration is 
the product of the original metric tensor and the cr- function e"^"", the equations given by 
the variations over g^^^, are the combinations of the equations given by the variation over 
the original metric and cx-equation. 

It often happens that we can drop the terms linear in a in (|^). In particular, one can 
redefine the corresponding source term as it is in the case of IR sector of 4D QG In 
the following, we only consider this case. Then the variations of 5* + Tind + r[l, gj^J] with 



respect to g 




i?4 




p and cr are given by (see also ^) 
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Ul 2 

idraf - dlo + 2dradrP\ + f-8e'^9,a9, fe-^^S^ 

J 1 n L V 



16 



32 



- W + + ^')| hQe^PdradrR^ + A{drafdlp 
I 6 J 16 L 



X 



]ly (-adla + {draf 



h" + \{h + h')^[dla + {draf) 



12 



-he^^dr \ In 



/?4 

7, 



dj-R^ 
3 



-[^d'^p-2{drpr-\d', + ldrpdr 



=-bd-^.p\n^ 



+ 



3 



p 



4 



3^ l^rPJ + 



d?.RA 



+ 



4 A 



+ b' dip 



1 

-R4 



-2e-2''92p + 
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-9,2 e'" + ^e 



2p+2a 



+ h'rl{-2{dlafe- 



2p 



-i {-2 (9.cr)2e-2^9,V + {{d^afe-'P)] 
■ {fe" + ^(6 + h')yi [dl [e-'P (dla + {d^af)] 



(9) 
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+^0 



In 



i?4 



4 /6 



+ b' dl In 



r. 



+ lte-2''9,2p9,2 jln 



\ 



-R4 
7, 



4 



69? In 



2 266^'' 1, 

3 



i?4 

4 /6 



-6e-2p [dip) + ^ + Tj65,2i?4 - I o + 1 ^'p 



' 



'0 



' 2q-^p { 2 



Ra 



be-^P{dlpf + 



2be'^p 1 



3"' ^"'^"^ ■ 3rg '3 



+-bdtRi- , 



4 /6 



3r5 6 



6" + ^(6 + 6')^o' 



^9^-R4 + ^ {9r-R49rCr + drRidrfj} 



+?>dl [e~^P (dla + {draf) } - 69, {e-^^S.a (dla + } 
Note that now the real 4 dimensional metric is given by 

ds" = -e^'^+^Pdx+dx- + r^dn' . 
Let us make the following change: dr etc. as 

dr = eP+'^di , dl = e^P^^'' {df + dipdi + diadi) . 
Using the new radial coordinate I, the metric (|12D is rewritten as follows: 



ds' 

f{l) 
r{l) 



-f{i)dt' + di' + r{iydn' 

roe" . 



(10) 



(11) 



(12) 



(13) 



(14) 



Here /(/) is called a redshift function and r(/) is a shape function. If /(/) and r{l) are 
smooth positive-definite functions, which satisfy the conditions: 



/(/) ^ 1, r(/) ^ / when |/| ^ 00 
/(O, r{l) finite when |/| ^ , 



(15) 



the metric expresses the wormhole which connects two asymptotically fiat universes. 
For the choice of the metric (p^, Eqs.(^, (^, (p!OD and (0) are rewritten as follows: 
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g4p+4o- 



-dfa + diadip 



16 



-Se-^'^diadi (e^- (af (7 + dipdia + {diaf)) 



3 

^2g4p+4o- 



16 



lQ>e-^''diadiR^ + 4(5^(7)' (Sfp + {dipf + S^aS^p) 

-12e-2'^9;aai (e^- (d^a + ^^p^^a + 2 {diaf)) 
+12 (af (7 + 9ip9ia + 2 (9,(7)') (a^a)"' 

^2g2p+2cr MS 1 15 1 



xe 



16 

1p+2a 



6' (-a (dfa + 9,p9,(7 + {diaf) + (9,(7)') 



- + ^(6 + 6')} (5f ^ + + 2 (a,a)') 



+^0 



-e'^'^'- {i^fp - ^(9,p)' + la,aa,p - + ^a,pa, - \d,adi\ 
xe'^^^'^ (afp + {dipf + a,(7a,p) in (^^ j 

+ \(fl + dipdi^diadi) {lng)j 

+ 16 (d^p + ((?,p)' + a,aa,p)' + 

+ \(d^ + dipdi + diadi)R^ 



i?4 





-2e-2- (afp+ (a,p)' + diudip) + ^ , 



2„2p 
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-e''^ (af + dipdi + a,aa,) e''^ + ^e''^ 



'^0 



+6Vo'e'^+^'^ |-2 {d^a + 9,p9,a + (a,a)') 

(af p + (9,p)' + diudip) ((7 (af (7 + a,p9,(7 + {diuf)) 

+ \^'^" (df + dipdi + diadi) (ae''^ (sf a + dipdia + {diaf 
-\ {-2 {diaf {d^p + (a,p)' + diadip) 



(16) 



(17) 
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+e 



-2cr 



[d^ + dipdi + diadi) {e"^ {diaf)}} 



-3 (d1a^dipdxa^2{dxaf) 



^-\)(plp^{dxpf ^d,ad,p)''\i. 



i?4 



+ 



4 
"3 



In 



-R4 



+ 



6') (af + 9,p9^ + a,(7a;) In 



i?4 



+-6 + {dipf + (af + dipdi + jln (^-| j | 

-^fee-^'^ (af + dipdi + a^aa,) je^'^ (af + dipdi + ^^^a,) jin 

2 



i?4 

f-r^ + -i^h (df + dipdi + diadi) i?4 



bc''^ {dfp + {dipf + diadip) 



-2. (q2 



26 1 

+^ + fe^" {df + dipdi + diadi) R, 



2b 

Ae^'' (b 



rl V3 



-2e^'^ {3af(7 + Adipdia + 6 (a^a)' + dfp + (a^p)'} + 



„2(T 



+6Vo2e2^ [262-^ (af + dipdi + a^a^/) (c^-^ + dipdia + (^^a)^ 
+^ (e^'^ {dfa + a^pa^a + {diaf) (dfp + {dipf + diudip) 
W {df + dipdi + diadi) (ac2- {dfp + {dipf + diadip) p 

-— ^ (afcT + dipdia + (5^(7)^) - —diRidia 



6" + -{b + b')yy+^'' [-- (9f + dipdi + i?4 + diRA<y 

+3 (af + a,pa, + diad) {e^- (af (7 + dipdia + 2 (^^a)') } 



The equations ([T6|) and (p!8D contain the 6th order derivatives of p, df p and the equations 
(|I^), (|TB|) and ( ]TPD contain the 4th order derivatives of a. Therefore we need to impose 
the initial conditions including 5th order derivative of p and 3rd order derivative of a: 

dfp\i=o = dfp\i=o = dfp\i=o = dfp\ 



dip\i=o ^ 
dia\i=o 



U=o 



dfa 



dfa\i=o = 



(20) 



The equations (|1^) and (|T8|) contain the 5th order derivatives of a from dfR^ but dfa 
always appears in the form dfadp. Therefore these terms vanish at / = 0. This tells 
that we cannot impose the initial condition dfa = since we cannot solve the differential 
equations with respect to dfa at Z = 0. 

We should also note that all the equations (|TB|), ( ]T^ ) and (|T^ are not the dynamical 
equations of motion but the combinations with the constraint. The constraint can be 
found by putting / = under the initial conditions (|20|), when the equations (^) and 
(p!8|) have the following form: 








b + b' + -b" + Sb'ao ] a, 



(4) 




146e-2<xo^^2^(6) 



(21) 



3e-2'^« 
SGvr 



2fee- 



-4cro 



In 



W[b + b' + y -2b'a^ 



a, 



(4) 




+ 4 (2 6 + 3 b") a, 



(4) 



:22i 



(23) 



Equations (^ 



and (23) will be compatible if o"o satisfies to the following equation: 

2 \ 







-224 6 (2 6 + 3 b") G vr In 



+e2'"o (54 6 - 30 6' + 81 b" + 40 b' ao) . 

Then we can choose the initial conditions as (pOD and (0) with arbitrary p\ 
under these initial conditions, we can use any two of equations 



(24) 
Then 



^), ([TID and (|T9D as 
dynamical equations of motion. Equation (24) admits the real solution for ctq if 

2 



b (112 6+168 6") e"5 
and from (|24D we have the restriction 



27 (2 6+3 b") 

2067 G TT In 



p^ ro^^ 



< -5 6' 



(25) 



1 

ao>- 



54 6 + 816" 



From equations (|T6D, (|T^ 
of p and a near / = 

P{1) 
a{l) 



40 |6'| 

and initial conditions (12 



(26) 

we get the following behavior 



Po + 



(Jo - 



26 + 26' + 36" + m'ao 
80640 G TT ro^ 6 (26 + 36"; 



/6 + 0{l^) 



-2(70 



96G7rro2 (26 + 36' 



-f + oir 



(27) 
(28) 
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Then for initial conditions (^Dj) a{l) and p{l) are decreasing functions near / = 0. 

For the numerical calculation as an example we will consider A/" = 4 SU{N) super- YM 
theory where (see, for example, [p!0[] ) 

Such theory became popular recently in relation with AdS/CFT correspondence. We take 
it as typical example of GUT. Note also that as it was mentioned in ref. the explicit 
choice of b" does not influence the equations of motion. 



Then from (|27| - p8D we get 

P{1) 
a{l) 



Po + 



-CTq 



10080G7rro^6 



/^ + 0(/«), 



'2(70 



and from 



^) we obtain 

{prof > 2 - 
{prof = 2 - 

{prof < 2 - 



b>0, (To > 21/10 
6>0, (70 = 21/10 

5gl6/5 

< 6 < 



IUGtt ln(^ 
Let us choose for numeric investigation, that 



8 21 

5<"°<To- 



(30) 
(31) 



(32) 
(33) 

(34) 



G = l, fi^ = 2, ro^ = 2//i2 = 1. 

The results of numerical calculation are given by Figures 1 and 2 for redshift and shape 
functions where = 5. 




Hence we showed that for above choice of initial conditions the inducing of primordial 
wormholes at the early Universe is rather unrealistic. The throat radius of wormhole 
quickly shrinks to zero. 

Let us now consider another choice of initial conditions 

dip\i=o = dfp\i=o = dfp\i=o = 0, dfp\i=Q ^ 

dia\i=o = dfa\i=o = dfa\i=o = . (35) 
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From equations (|TB|), (]TB|), ( ]TP| ) and initial conditions ( P5| ) we get the following behavior 
of p and cr near / = 



1 



79262G 



TT 



-ggfee-^^" + 462G'7rpofXo (-51 + SGe^'^" + 20e^'"Vof^o^ 



4£v^Pof^o-\/fr'^G (gge^'^o (21 - IOctq) + feG'7r(-51 + See^-^o + 20e4'^opo^Q 

1 



(36) 



132 62 G 



6^ G (51 - 36 e^"° - 20 e^"" po a^' 



+ eJb^G (gge^'^o (21 - lOao) + ^^vr (-51 + 36e4'^o + 20e4'^o 



(37) 



where e = ±1. 

Moreover we obtain the conditions 



dfp\i=o = d^cT\i=o = 0, 



or 



Sgie^-^" (10 (To -21) = SfeCvr (-51 + 36e^'"o + 20e^'"Vof^o^)', (^o > 



21 
10' 



(3^ 



(39) 



For case ao = 21/10 we have 

1 /-275e-2i/5 



dfa\i=o 
dfp\i=o 



2200 
1 



110 



V bGn 
85 - 60e^2/5-147e^2/5^^ 



+ (1680 e^2/5 _ 2380) po + 4116 e^^/s 



(40) 
(41) 



Then functions a{l) and p(/) will increase near / = if 

-42/5 



Po < 



2g4 6G' 



bG (-12 + 17e-^'/^) - 



bG (231ef + b(l7 -Ueff Gn^ 



for (7 = 1 and b = {N'^ — l)/(87r)^ we have po < —0.4081. The explicit numerical study 
(see Figures 3 and 4) shows that shape and redshift functions are increasing giving the 
window for inducing of primordial wormholes. 



8.4 



8.3 



8.2 



Fig. 3 n=10 



Fig. 4 n=10 



0.2 0.4 0.6 0.8 

On these figures (Tq = 21/10 and po = —0.41. 




0.2 0.4 0.6 
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Hence we numerically proved that at least for some initial conditions GUTs at the early 
Universe may help in producing of primordial wormholes. It is of course the open question 
which initial conditions from the two classes presented above lead to more stable config- 
uration. Moreover, due to complicated structure of field equations and initial conditions 
itself we cannot classify from the very beginning the initial conditions as supporting (or 
not) wormholes production. Nevertheless for any specific GUT under discussion above 
study may be easily repeated and principal possibility of wormholes inducing may be 
shown at least numerically. 
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